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Abstract 

Electromagnetic  material  characterization  is  the  process  of  determining  the 
complex  permittivity  and  permeability  of  a  material.  A  new  waveguide  material 
measurement  technique  is  developed  to  reduce  test  sample  size  requirements  for  low- 
frequency  applications.  Specifically,  a  waveguide  sample  holder  having  a  reduced 
aperture  is  utilized  to  decrease  the  time  and  cost  spent  producing  large  precision  test 
samples.  This  type  of  sample  holder  causes  a  disruption  in  the  waveguide-wall  surface 
currents  that  results  in  the  excitation  of  higher-order  modes.  This  thesis  will  demonstrate 
how  these  higher-order  modes  can  be  accommodated  using  a  modal-analysis  technique, 
thus  resulting  in  the  ability  to  measure  smaller  samples  mounted  in  large  waveguides  and 
still  determine  the  constitutive  parameters  of  the  materials  at  the  desired  frequencies. 

Experimental  results  for  acrylic  and  magnetic  and  electric  radar  absorbing 
materials  are  compared  to  ideal  full-aperture  measurements  at  S-band  (2.6-3.95  GHz) 
to  verily  the  modal-analysis  technique.  In  addition,  since  uncertainties  in  test  sample 
thickness  and  placement  are  predominantly  encountered  in  material  characterization 
measurements,  a  differential  error  analysis  is  performed  to  determine  associated  error 
bounds. 


IX 


STEPPED  WAVEGUIDE  ELECTROMAGNETIC  MATERIAL 


CHARACTERIZATION  TECHNIQUE 

1.  Introduction 

An  electromagnetic  material  is  characterized  by  its  permittivity  and  permeability. 
In  general,  the  permittivity  of  a  dielectric  material  determines  its  susceptibility  to 
becoming  polarized  when  exposed  to  an  electric  field.  Similarly,  the  permeability  of  a 
magnetic  material  describes  its  susceptibility  to  becoming  magnetized  when  exposed  to  a 
magnetic  field.  Permittivity  and  permeability  are  generally  both  complex  quantities, 
where  the  real  parts  are  related  to  energy  storage  and  the  imaginary  parts  are  related  to 
the  loss  mechanisms  of  the  material  that  convert  incident  electromagnetic  radiation  to 
heat. 

A  multitude  of  applications  require  measurement  of  the  constitutive  parameters  of 
an  electromagnetic  material.  For  example,  consider  the  following  overly  simplified 
process  description:  1)  electrical  engineers  specify  specific  constitutive  parameters  for  a 
new  radar  absorbing  material  (RAM)  coating  for  a  stealth  aircraft,  2)  chemical  engineers 
labor  to  fabricate  the  RAM,  3)  to  test  the  material  formulation,  the  responsibility  returns 
to  the  electrical  engineers,  and  4)  the  process  continues  until  an  acceptable  formulation  is 
attained  or  the  original  specification  is  modified.  Meanwhile,  during  this  often  very  time- 
consuming  process,  large  costs  are  incurred.  This  thesis  offers  a  waveguide  material 
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measurement  teehnique  aimed  at  redueing  the  test  sample  size  requirements  for  low- 
frequency  applications,  and  should  thus  lead  to  significant  time  and  cost  savings. 

1.1  Problem  Statement 

Typically,  the  material  characterization  process  involves  several  basic  steps. 

First,  a  test  sample  is  carefully  machined  to  fit  into  the  testing  device.  Next,  the  device  is 
connected  to  a  vector  network  analyzer  to  obtain  the  experimentally  measured  scattering 
parameters.  The  network  analyzer  operates  by  launching  an  incident  wave  towards  the 
test  sample  and  subsequently  measuring  the  amount  of  signal  that  is  reflected  from  and 
transmitted  through  the  material.  Finally,  numerical  algorithms  are  used  to  extract  the 
complex  permittivity  and  permeability  of  the  test  sample  from  the  experimental 
S-parameters.  The  predominant  algorithm  used  to  compute  the  material  parameters  is  the 
Nicolson-Ross-Weir  (NRW)  technique  [l]-[5]. 

The  NRW  technique  has  the  advantage  of  being  one  of  the  few  “direct”  extraction 
methods  where  the  material  parameters  are  calculated  analytically  in  closed-form  from 
the  experimental  S-parameters.  Most  other  extraction  methods  are  “iterative”  in  nature 
with  an  initial  guess  required  to  iteratively  find  the  material  parameter  solutions  that 
satisfy  an  optimization  algorithm. 

Test  samples  must  be  linear,  homogeneous,  and  isotropic  with  coplanar  front  and 
back  surfaces  for  proper  application  of  the  NRW  technique.  In  addition,  only  a  single 
mode  is  assumed  to  be  present  in  the  material.  Thus,  for  waveguide  applications,  the 
material  must  completely  fill  the  waveguide  cross-section  to  prevent  the  excitation  of 
higher-order  modes.  For  frequencies  greater  than  several  gigahertz,  producing  test 
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samples  is  relatively  simple.  In  an  X-band  (8.2-12.4  GHz  )  waveguide,  for  example, 


test  sample  dimensions  in  the  eross-sectional  plane  are  only  0.9  by  0.4  inehes.  However, 
for  lower-frequency  applications  waveguide  dimensions  become  progressively  larger. 
Consequently,  larger  quantities  of  materials  are  required  leading  to  possible  sample 
fabrication  difficulties. 

To  overcome  these  large  test  sample  fabrication  difficulties  imposed  by  the  NRW 
technique  for  low-frequency  waveguide  measurements,  solutions  in  the  past  have  been  to 
only  partially  fill  the  cross-section  of  the  sample  holder  and  perform  an  analysis  using 
higher-order  longitudinal  section  electric  (LSE)  or  longitudinal  section  magnetic  (LSM) 
propagating  modes  [6], [7].  While  the  LSE/LSM  approach  reduces  the  test  sample  size 
requirements,  it  suffers  three  main  drawbacks.  Eirst,  the  higher-order  mode  propagation 
constants  must  be  found  numerically  by  solving  transcendental  equations.  Second,  it  is 
difficult  to  associate  the  propagation  constants  with  the  correct  modes.  And  third,  this 
iterative  method  is  very  sensitive  to  the  initial  guesses  for  the  constitutive  parameters. 

The  goal  of  this  thesis  is  to  eliminate  the  difficulties  in  dealing  with  a  partially 
filled  waveguide  cross-section  by  completely  fdling  a  sample  holder  having  a  reduced 
aperture.  Unlike  the  ESE/ESM  approach,  it  will  be  shown  that  the  reduced-aperture 
sample  holder  leads  to  well-known,  closed-form  expressions  for  the  propagation 
constants  and  their  associated  modes.  Specifically,  the  higher-order  modes  that  result 
from  the  disruption  of  the  waveguide-wall  surface  currents  when  encountering  the 
waveguide-feed/sample-holder  junctions  will  be  accommodated  through  a  modal-analysis 
technique. 
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1.2  Limitations 


While  the  modal-analysis  technique  is  much  easier  to  implement  than  the 
LSE/LSM  approach  due  to  its  closed-form  solutions,  there  are  limitations  that  must  be 
understood.  Specifically,  high-loss  materials  measured  in  greatly  reduced  apertures  will 
lead  to  transmission  coefficients  that  approach  the  noise  floor  of  the  network  analyzer. 
The  poor  transmission  coupling  will  in  turn  significantly  degrade  the  material  parameter 
extraction  process.  Thus,  with  high-loss  materials  of  greatly  reduced  test  sample  sizes, 
partially  filling  a  full-aperture  sample  holder  and  accounting  for  the  higher-order 
LSE/LSM  modes  will  lead  to  better  transmission  responses,  but  at  the  cost  of  solving 
complicated  transcendental  equations.  As  with  any  engineering  problem,  there  are 
tradeoffs  to  using  different  approaches. 

1.3  Scope 

Many  other  measurement  devices  can  be  utilized  for  material  characterization 
(including,  but  certainly  not  limited  to,  coaxial  transmission  lines,  free-field  measurement 
systems,  strip  transmission  lines,  microstrip  field  applicators,  and  cavity  resonators),  but 
this  thesis  will  explore  only  rectangular  waveguide  applications.  In  addition,  only  S-band 
(2.6-3.95  GHz)  measurement  results  will  be  presented  as  a  proof-of-concept 
demonstration.  However,  the  stepped  waveguide  modal-analysis  technique  will  be 
developed  in  a  general  manner  that  is  fully  adaptable  to  the  low-frequency  measurements 
for  which  it  was  designed. 
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1.4  Thesis  Organization 


Chapter  2  provides  a  theoretical  background  focusing  on  transverse  electric  (TE) 
modes  in  rectangular  waveguides  and  wave  transmission  matrices  for  analysis  of  multi¬ 
layer  environments.  Chapter  3  derives  the  modal-analysis  technique  that  allows  a  smaller 
test  sample  to  be  mounted  in  a  large  waveguide  and  still  determine  the  constitutive 
parameters  of  the  material  at  the  desired  frequencies.  Chapter  4  presents  several 
experimental  results  that  verify  the  modal-analysis  technique.  In  addition,  a  differential 
error  analysis  related  to  uncertainties  in  measuring  test  sample  thickness  and  placement 
within  the  sample  holder  is  conducted  in  Chapter  4  to  form  error  bounds  on  the  extracted 
material  parameters.  Finally,  a  brief  set  of  conclusions  and  recommendations  for  future 
research  are  given  in  Chapter  5. 
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2.  Background 


This  chapter  provides  the  baekground  neeessary  for  understanding  the  modal 
analysis  method  used  in  this  thesis  to  aecommodate  the  stepped-waveguide  structure. 
Specifically,  modes  of  rectangular  waveguides  and  wave  transmission  matrices  will  be 
briefly  reviewed.  Although  some  minor  speeializations  have  been  made,  an  experienced 
reader  familiar  with  both  of  these  eoneepts  should  be  able  to  skip  to  the  next  ehapter 
without  loss  of  understanding. 

2.1  TE  Modes  of  Rectangular  Waveguides 

Perfectly  conducting  rectangular  waveguides  can  support  Transverse  Eleetric 
(TE)  and  Transverse  Magnetic  (TM)  modes.  Most  praetieal  waveguide  systems  are 
designed  to  operate  so  that  only  the  dominant  TEio  propagating  mode  is  exeited  because 
it  is  the  lowest-order  mode  and  thus  prevents  the  possibility  of  dispersion  from  higher- 
order  propagating  modes.  Since  TEio  mode  excitation  was  used  for  all  laboratory 
measurements  in  this  thesis,  discussion  will  be  limited  to  developing  the  pertinent  field 
expressions  for  TE  modes.  The  interested  reader  ean  obtain  a  more  general  review  of 
guided-wave  theory  in  Collin  [6]. 

2,1.1  Helmholtz  Wave  Equation. 

The  most  fundamental  beginning  to  any  electromagnetic  analysis  is  Maxwell’s 
Equations.  Eor  a  linear,  homogenous,  isotropic  material  in  a  souree  free  region,  the  time 
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harmonic  representation  (with  a  suppressed  time  dependenee)  of  Maxwell’s 
equations  are  given  by 


V  X  E  =  -jeopH 

(2.1) 

V  X  //  =  jeosE 

(2.2) 

V-E  =  0 

(2.3) 

V//  =  0 

(2.4) 

where  E,  H  are  the  eleetrie  and  magnetie  field  intensity  veetors  whieh  are  funetions  of 
position  and  angular  frequency  co ,  and  6’,//  are  the  effective  permittivity  and 
permeability  of  the  medium  which  are,  in  general,  also  eomplex  functions  of  angular 
frequency  co .  To  uncouple  Maxwell’s  equations,  the  vector  identity 

VxVxK  =  V(V-  V)-V^V  (where  V  is  an  arbitrary  vector)  is  used  to  obtain  the  time 
harmonie,  homogenous  Helmholtz  vector  wave  equations 

+  1^1  =  0  (2.5) 

where  the  complex  wave  number  k  =  co^[^ . 

In  dealing  with  wave-guiding  struetures,  it  is  eonvenient  to  break  up  the 
Helmholtz  equation  into  its  transverse  and  longitudinal  eomponents  for  ease  of 
identifying  TE  and  TM  mode  solutions.  To  start,  the  eleetrie  and  magnetic  fields  must 
also  be  separated  into  their  respeetive  transverse  and  longitudinal  parts.  For  forward  and 
reverse  propagating  waves  along  the  z-axis,  the  fields  can  be  expressed  as 

E{r)  =  E{p,z)  =  E^{p,z)±zE^{p,z)  =  e^(p)e+^^  ±%(p)e+^^  (2.6) 
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(2.7) 


H  (r  )  =  H{p,  z)  =  ±H^  {p,  z)  +  zH^  {p,  z)  =  ±h^  (p)e^^"  +  zh^  (p)e^^" 
where  7  is  the  propagation  constant,  the  position  vectors  r  =  (x,  7,  z)  and  p  =  (x,  7) , 
and  the  subscripts  t,  z  are  used  to  represent  the  transverse  and  longitudinal  components, 
respectively  [8].  For  TE  modes,  in  which  e^=0,  Maxwell’s  equations  reduce  to  the 
following  expressions  upon  equating  the  transverse  and  longitudinal  field  components 


X 

II 

1 

(2.8) 

V^xe^  =  -zjcoph^ 

(2.9) 

+  /zxhf  =  -jcose^ 

(2.10) 

0 

II 

X 

(2.11) 

II 

0 

(2.12) 

^t-ht=rK 

(2.13) 

where  V  =  +  z  and  the  vector  identity  VxsV  =  VsxV  +  sV  xV  (where  s  is  an 

arbitrary  scalar)  are  used.  Using  the  vector  identity  AxBxC  =  b(^A-C^- A- ,  the 

transverse  field  components  can  also  be  expressed  in  terms  of  the  longitudinal 
component  as 


h,=-T2^A 

(2.14) 

et=-ZTE^X^ 

(2.15) 

where  the  TE  wave  impedance  is  given  by 
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Zte  -■ 


(2.16) 


Using  ,  the  wave  equation  for  H  can  now  be  decomposed  into  its 


transverse  and  longitudinal  components 


W^h,{p)  +  klh,{p)  =  Q 


(2.17) 


V,\(p)  +  ^,\(/5)  =  0 


(2.18) 


where 


(2.19) 


2,1.2  Separation-of-Variables  and  Boundary  Condition  Enforcement. 

It  has  just  been  shown  that  if  is  known,  all  the  field  components  of  TE  modes 
can  be  generated  using  (2.14)  and  (2.15).  Therefore,  solution  to  the  scalar  wave  equation 
for  given  by  (2.18)  will  be  sought  using  the  separation-of-variables  technique.  With 

the  product  solution  h^i^'p)  =  assumed,  (2.18)  becomes 


Y(y)  +  X{x)  +  k^X(x}Y(y}  =  0 

dx  dy 


(2.20) 


and  then  dividing  by  X(x)Y(y)  gives. 


1  d^Xjx)  1  d^Yjy)  2 

X(x)  dx^  Y(y)  dy^ 


(2.21) 


where  the  definitions  k^,ky  are  made  observing  that  each  term  of  (2.21)  must  be  equal  to 
a  constant,  and  it  is  evident  that  k^  =k^  +  k^  .  The  solutions  to  the  two  separate  ordinary 


differential  equations. 
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(2.22) 


d^X{x)  ,2^,  ,  ^ 

dx^ 

dy 

are  found  considering  the  rectangular  waveguide  geometry  presented  in  Figure  2-1,  with 
the  origin  shifted  to  the  cross-sectional  center  for  mathematical  convenience  in  the  next 
chapter.  With  the  wave  bounded  in  x  and  y  ,  the  appropriate  solutions  to  (2.22)  imply 
the  following  convenient  choice 


Kixy)  =  X{x)Y{y)^  sin { ( X  +  ^)}  +  ^2  cos { ( X  +  ^ )} 


Cl  sin \ky  [y  +  %)}  +  C2  cos  (>^  +  %)} 


.  (2.23) 


Imposing  the  well-known  boundary  condition  that  the  tangential  component  of  the 
electric  field  at  a  perfectly  conducting  surface  must  be  zero,  expressed  mathematically  as. 


nxE  =  0 


hxe=0  ...  since  e=0 


(2.24) 


it  can  be  shown  with  the  help  of  (2.14)  and  (2.15)  that,  at  the  waveguide  walls, 

.  dhyix.y) 

n-V^h^{x,y)  = - ^ - 2  =  o 

on 


(2.25) 


where  h  is  the  unit  normal  vector  directed  into  the  field  region  from  the  waveguide 
walls.  Enforcing  the  boundary  conditions  implied  by  (2.25),  the  generating  function  for 
TE  wave  propagation  in  a  rectangular  waveguide  is  given  by 


K  ix  y)  =  cos  {k^  (x  +  ^)}  cos  [k^,  (y  +  %)} 


k^  = - ,  m  =  0,1,2... 

a 

k  =—,  n  =  0,1,2... 
b 


(2.26) 
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A  y 


A'  z 

Figure  2-1 .  Geometry  of  a  reetangular  waveguide  with  the  origin  plaeed  in  the  eross- 
sectional  eenter. 


where  is  the  mn*  mode  amplitude.  It  must  also  be  noted  that  m  and  n  eannot  both 

be  zero  simultaneously,  otherwise,  upon  examination  of  (2.8)-(2.15),  a  trivial  solution  of 
no  wave  is  reaehed. 

Finally,  applying  the  generating  funetion  expressed  in  (2.26)  to  (2.14)  and  (2.15) 
gives  the  transverse  field  equations 


G  =  ^  AAy  COS (x  +  ^)} sin [y  +  Yi)] 
G  =  sin  {k,  [x  +  ^)}  cos  [k^,  [y  +  YY)] 

K=-— 

Z 

^TE 

hy  =  — 

^TE 


(2.27) 


2,1.3  Cutoff  Frequency. 

The  cutoff  frequency  of  the  mn*  rectangular  waveguide  mode  is  defined  by 
setting  the  propagation  constant  of  (2.19)  to  zero: 
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Ymn  =  =  0 


(2.28) 


Thus,  rearrangement  of  (2.28)  gives  the  eutoff  frequency  of  the  TEmn  mode  as 


(  fco  }  ffifi 


+ 


mn 
a  j 


nn 

K  ~b  j 


m  =  0,1,2,...] 

n  =  0,l,2,...J 


m  =  n^0 


(2.29) 


For  example,  consider  a  free-space  fdled  S-band  waveguide  with  dimensions 

a  =  2.84  inches  and  b  =  \34  inches .  The  cutoff  frequency  for  the  dominant  TEio  mode 

is  2.078  GHz  ,  while  that  of  the  next  higher-order  mode,  the  TE20  mode,  is  4.156  GHz . 

S-band  waveguides  are  customarily  operated  in  the  frequency  range 

1.25  10  <  /  <  0.95 20  2.6  </  <  3.95  GHz  ).  This  operational  regime  ensures  a 


well-behaved  dominate  propagating  mode  with  all  higher-order  modes  being  evanescent. 


2.2  Wave  Transmission  Matrices  (A-Parameters) 

Wave  transmission  matrices  offer  a  simple  method  to  deal  with  multi-layered 
environments,  and  will  be  utilized  in  the  next  chapter  to  model  the  sample  holder  region. 
This  section  will  develop  individual  2x2  matrices  that  describe  the  relationship  between 
incident  and  reflected  wave  amplitudes  at  a  prescribed  output  terminal  plane  to  those  at 
the  input  terminal  plane.  These  individual  matrices  are  then  easily  cascaded  to  obtain  the 
overall  wave  amplitudes  relationships  of  a  multi-layered  environment  [6]. 

2,2.1  Geometry  of  Multi-Layered  Environment. 

Consider  the  multi-layered  system  environment  of  Figure  2-2  comprised  of  N 
material  layers.  The  respective  complex  permittivity,  permeability,  and  thickness  of  the 
material  layer  are  £•. ,  ju. ,  and  f .  The  complex  wave  amplitudes  of  the  incident  and 
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reflected  waves  within  the  layer  are  represented  by  ,  respectively.  Finally,  r.,r; 
correspond  to  the  Fresnel  interfacial  reflection  and  transmission  coefficients  at  the  input 
plane  of  the  layer. 

2.2,2  A-Parameter  Description  of  Single  Layer. 

To  analyze  the  multi-layered  system  presented  in  the  previous  section,  it  is 
necessary  to  first  examine  the  expanded  view  of  a  single  layer  as  shown  in  Figure  2-3.  If 
r2,r2and  are  the  respective  interfacial  reflection  and  transmission  coefficients 


experienced  by  waves  Cj  and  b'^ ,  then  the  following  relations  prevail 


h 

c'l 


FjCj  +  T2b2 
^  ^2^2  ^2^1 


(2.30) 


Then,  using  the  following  well-known  Fresnel  equations  (based  on  continuity  of 
tangential  electric  and  magnetic  fields  at  the  interface) 

^2  “  1  ^2 

r;  =  i+r;  (2.3 1) 

7  -7 

r2  =  -r;  =  ^ ^ 

Z2  +  Z, 


and  a  little  algebra  leads  to  the  matrix  expression 


Cl 

1 

■  1 

r2' 

C2 

A  J 

^2 

^2 

1 

K_ 

(2.32) 


where  Z.  is  the  wave  impedance  of  the  layer.  Thus,  (2.32)  describes  the  relationship 
between  the  forward  and  reverse  traveling  waves  (cj,hj)  immediately  to  the  left  of  the 
interface  to  the  forward  and  reverse  traveling  waves  immediately  to  the  right  of 

the  interface. 


2-8 


^1 

<-2 

1 

A^-1 

^2 

w 

1 

1 

1  T 

i 

^  W 

^W-l 

^  W 

i 

1 

i 

Figure  2-2.  Multi-layer  environment. 


Figure  2-3.  Expanded  view  of  a  single  layer  used  to  demonstrate  the  refleeted  and 
transmitted  waves  at  a  planar  interface  as  well  as  the  propagation  phase  delay. 


To  complete  the  matrix  description  of  the  layer  under  consideration,  the  waves 
{€'2,^2)  must  be  related  to  waves  {c2,b2)  after  propagating  a  distance  Assuming  a 

linear,  homogeneous,  and  isotropic  medium,  the  complex  wave  amplitude  relationships 
are 


c,  =  c\e 

c;l 

\ 

0 

1 

<N 

1 _ 

1 

II 

b\\ 

0 

— 1 

1 

(2.33) 
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where  Yi  is  the  propagation  constant  in  layer  2.  Finally,  combining  the  results  of  (2.32) 


and  (2.33)  yields 


Cl 

1 

erih 

C2  ! 

“Al 

1 - 

<N 

^2 

a: 

u 

e-nfi 

1 — 

ts> 

1  _ 

_ 1 

^22  _ 

A 

(2.34) 


where  the  transfer  function  is  normally  described  in  terms  of  A-parameters . 

2.2,3  A-Parameter  Description  of  Multi-Layered  Environment, 

Using  the  single-layer  A-parameter  description  developed  in  the  previous  section, 
a  generalized  result  can  be  applied  to  the  multi-layered  system  of  Figure  2-2.  The  overall 
A-parameter  description  of  an  A-layer  system  is  given  by  the  cascaded  product  of  the 
wave  matrices  of  each  individual  layer 


N  1 

=nf 

i=2 


JiZi 


~Yi^i 


•'N 


4 1 

A  A 
A21  ^22. 


(2.35) 


where 
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r,=^ — ^  ,  T=i+r, 


(2.36) 


2.3  Summary 

This  chapter  reviewed  TE  modes  of  rectangular  waveguides  and  wave 
transmission  matrices.  Understanding  of  both  concepts  will  form  the  basis  of  the  modal- 
analysis  technique  presented  in  the  next  chapter.  It  will  be  shown  that  under  TEio  mode 
excitation,  the  geometry  chosen  for  the  reduced-aperture  sample  holder  will  cause  only 
specific  TE  modes  to  be  scattered  -  no  TM  modes  will  need  to  be  considered.  Also,  the 
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wave  transmission  matrix  will  provide  a  eonvenient  method  to  handle  measurements  of 
multiple  material  layers  within  the  sample  holder  region. 
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3.  Methodology 


This  chapter  presents  a  modal  analysis  proeess  that  can  be  used  to  extraet  the 
constitutive  material  parameters  of  a  test  sample  mounted  in  a  stepped-waveguide 
structure  from  the  experimentally  measured  S-parameters,  5*“^ ,  of  a  network  analyzer. 
The  modal  analysis  proeess  will  be  broken  up  into  three  simple  steps:  1)  expand  the  fields 
for  each  region,  2)  satisfy  boundary  conditions  across  each  interface,  and  3)  test  with 
appropriate  integral  operators. 

The  outputs  of  the  modal  analysis  process  are  the  theoretieally  calculated 
S-parameters,  Sf^ ,  for  an  estimate  of  the  eomplex  relative  permittivity  and  permeability 

of  the  test  sample.  Using  an  iterative  nonlinear  least-squares  optimization  algorithm,  the 
test  sample's  material  parameters  will  be  those  that  match,  within  a  specified  accuracy, 
the  theoretical  S-parameters  to  the  experimental  values  at  eaeh  frequency: 


min 


s,/^) 


(3.1) 


where  /  =  1 1, 2 1, 22, 12  . 


3.1  Field  Expansions 

The  stepped-waveguide  geometry  was  carefully  chosen  to  simplify  the  mode¬ 
matching  process  aeross  the  waveguide  to  sample  holder  interfaees.  By  matching  the 
reduced  aperture  of  the  sample  holder  to  the  shape  of  the  TEio  mode  excitation,  mode 
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Figure  3-1.  Top  view  stepped-waveguide  fixture  showing  the  sample  holder  region  as  an 
N-layer  environment. 


orthogonality  was  used  to  greatly  decrease  the  complexity  of  the  field  expansions  for  the 
different  scattered  modes. 

3,1,1  Stepped-Waveguide  Geometry, 

Consider  the  geometries  of  Figures  3-1  and  3-2  showing  respective  top  and  side 
views  of  the  stepped-waveguide  structure.  The  origin  is  placed  in  the  cross-sectional 
center  of  the  waveguide  for  mathematical  convenience.  Regions  A  (z  <  O)  and  B 

(z  >  are  air-filled  waveguide  sections  having  width  a  and  height  b  .  The  reduced- 
aperture  sample  holder  attached  between  regions  A  and  B  consists,  in  general,  of  a  multi- 
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layered  material  stack.  The  sample  holder  is  assumed  to  have  a  reduced  width  t 
symmetrically  centered  on  the  origin,  but  no  reduction  in  height.  The  thickness  of  each 
layer  in  the  sample  region  must  comprise  the  total  sample  holder  thickness 

N 

(3.2) 

1=1 

3.1,2  Mode  Orthogonality. 

The  disruption  in  the  waveguide-wall  surface  currents  experienced  at  the 
waveguide  to  sample  holder  junctions  will  cause  the  excitation  of  higher-order  modes. 
Examining  (2-27),  the  transverse  fields  of  the  TEio  mode  are  readily  seen  to  be 
y-invariant  and  even  about  the  center  width  of  the  waveguide  (x  =  O) .  Thus,  since  the 

aperture  of  the  sample  holder  is  reduced  only  in  width  and  centered  about  x  =  0  ,  all 
higher-order  excited  modes  will  also  be  y-invariant  and  even  about  the  waveguide  width 

[9]. 

Considering  only  this  limited  set  of  possible  higher-order  mode  shapes,  the  types 
of  higher-order  excited  modes  is  significantly  decreased.  Although  not  specifically 
shown  here,  all  TM  modes  are  y-dependent  and  can  thus  be  immediately  discounted  [8]. 
In  addition,  y-invariance  of  TE  modes  implies 

=0 

n  =  0  ^  k,,  =0  ^  ^  (3.3) 

>■ 

which  permits  only  TEmo  modes  with  the  field  equations 

(3.4) 
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Then,  using  a  trigonometric  identity,  the  following  simplification  can  be  made 


sin{^,(x  +  ^)} 


(-1)  cosk^x 
(-1)’^  sinA:^x 


,  m  =  l,3,5,... 

,  m  =  2,4,6,... 


(3.5) 


Finally,  the  even  requirement  about  x  =  0  leaves  only 


JCDjU 


A^„k^  cos  kx  =  cos  kx 


2  mn  X 


h=- 


y 

z 


m  =  1,3,5, 


(3.6) 


where  has  been  introduced  to  consolidate  the  constants,  and  the  subscript  on  the 

wave  impedance  has  been  suppressed  since  only  TE  modes  are  possible. 

Note  that  the  physical  insight  leading  to  consideration  of  only  this  limited  set  of 
higher-order  modes  can  be  rigorously  proven  through  mode  orthogonality.  That  is, 
integrating  the  dot  product  of  the  transverse  fields  of  the  TEio  mode  with  any  other  mode 
over  the  cross-section  of  the  sample  holder  aperture  will  yield  zero.  Eor  example, 
consider  integration  over  the  sample  holder  aperture  of  the  dot  product  of  the  electric 
fields  of  the  TEio  mode  with  the  TEmo  modes  for  even  values  of  m: 

bjl  tjl 

J  J  ^TEia'^TE  ’  m=  2,4,6...  (3.7) 

-6/2  -i/2 

Dropping  the  constants  and  ignoring  integration  over  y  (since  all  TEmo  modes  are 
y-invariant)  leads  to 

nx  .  mnx  ,  ^  ^ 

cos  —  sm - ax  =  0  ,  m  =  2,4,6...  (3.8) 

a  a 


1 
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which  is  readily  identified  as  the  integration  of  an  odd  function  over  symmetric  limits 
yielding  zero.  Similarly,  mode  orthogonality  can  be  used  to  disregard  all  TEmn  (n  ^  0) 
and  TM  modes.  Thus,  the  only  modes  that  will  couple  into  the  sample  holder  region  will 
be  TEmo  modes  for  odd  values  of  m. 

3.1,3  Field  Expansions  in  Waveguide  Regions, 

The  transverse  fields  in  the  waveguide  regions  are  given  by  the  following 
expansions 


Ef  =  ate,e  +  V  a  e  e^"" 

t  1  1  m  m 

m=l 

_  _  M  ^ 

Hf  =  -V  a.  h  e^ 

t  1  1  m  m 


...  z  < 0 


(3.9) 


Ef  =^bl,e^e 


m=\ 

N 


n=\ 


...  z>  d. 


(3.10) 


where  represents  the  eomplex  amplitude  of  the  incident  TEio  mode  traveling  in  the 
positive  z-direetion,  represents  the  eomplex  amplitude  of  the  TEmo  mode  reflected 


from  the  waveguide  to  sample  holder  junction  at  z  =  0 ,  hj,  represents  the  complex 
amplitude  of  the  TEmo  mode  transmitted  through  the  sample  region  at  z  =  ,  and  the 

negative  sign  in  the  Hf  equation  represents  the  polarization  rotation  upon  refieetion, 

which  is  consistent  with  the  direction  of  power  transport  for  the  refieeted  wave.  In 
addition,  the  following  indexing  eonvention  has  been  adopted 
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m  =  \,2,...,M  (3.11) 


e„=ycosk^^x  ,  k^^={2m-l)^y 


= 


zxe„ 


= 


/m 


a 

/m  = 


-kl 


where  k^=o)lc  with  c  =  l/ ^JsJU^  as  the  speed  of  light  in  a  vacuum. 

While  an  infinite  number  of  higher-order  modes  are  actually  excited  at  the 
waveguide  to  sample  holder  junction,  the  total  number  of  modes  M  that  must  be 
considered  for  a  given  convergence  threshold  is  dependent  on  the  degree  of  reduction  to 
the  sample  holder's  aperture.  Analogous  to  Fourier  Transform  theory,  where  a  narrow 
pulse  in  the  time  domain  translates  to  high-frequency  components  in  the  frequency 
domain,  larger  disruptions  to  the  waveguide-wall  surface  currents  will  require  more 
higher-order  modes  to  accurately  represent  the  discontinuity.  For  example,  in  an  ideal 
sample  holder  that  matches  the  dimensions  of  the  waveguide,  only  the  dominant  TEio 
mode  need  be  considered.  However,  as  the  width  reduction  increases,  there  is  an 
appreciable  amount  of  energy  shifted  to  successive  higher-order  modes.  Thus,  as  the 
sample  holder  aperture  decreases,  more  higher-order  modes  are  necessary  to  extract  the 
correct  material  parameters. 

3,1.4  Field  Expansions  in  Sample  Region, 

Modeling  the  sample  holder  as  a  reduced-aperture  multi-layer  environment  has 
two  important  advantages.  First,  it  allows  for  an  extra  degree-of- freedom  in  the 
manufacturing  process  by  not  requiring  the  test  sample  to  completely  fill  the  sample 
holder  thickness.  In  addition,  for  very  low  frequency  applications  where  waveguide 
dimensions  become  very  large,  a  rigid  dielectric  sandwich  can  be  used  to  support  thin, 
flimsy  test  samples.  However,  it  is  assumed  that  only  one  layer  in  the  sample  region  has 
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unknown  material  properties.  Thus,  in  the  case  of  a  dielectric  sandwich,  each  dielectric 
support  would  have  to  be  measured  separately  from  the  test  sample. 

The  transverse  field  equations  in  the  layer  of  the  sample  holder  region  are 
given  by 


M  M 


t  m,i  m 


+  75  e^e 


=1 

M 


n=\ 

M 


m=\ 


m=l 


\  d._,  <z<d, 


(3.12) 


where  the  superscript  5  is  used  to  indicate  the  sample  holder  region,  and  the  distance  to 
each  reference  plane  is  determined  by 

(3.13) 

k=\ 


An  indexing  convention  similar  to  (3.1 1)  has  also  been  adopted: 


e;=ycosO  , 


h!:'  = 


zxe 


= 


J(OMi 

yS,i 
/  m 


Ym  = 


‘  =  E: 


-kf 


m  = 


(3.14) 


where  the  complex  relative  permittivity  and  permeability  of  each  layer  determine  the 


wave  number  of  the  corresponding  medium  k.  • 


3.2  Satisfying  Boundary  Conditions 

The  second  step  of  the  modal  analysis  process  is  to  satisfy  boundary  conditions 
across  each  interface.  Enforcing  continuity  of  the  tangential  fields  across  the  waveguide 
to  sample  holder  interfaces  at  z  =  0,  leads  to 
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(3.15) 


M 


a,e,+  ?  i 


y  <  +  y  . . .  |x|  <  til 

m,[  m  w,l  m  II' 

m=\  m=l 

[  0  ...  -ajlKxK-tjl  ,tl2<x<al2 


iVi  ivi  ivi 

jI  hy-y  aji^  =  £  . . .  |x|  <  ^2  (3.16) 


(3.17) 


(3.18) 


ff2=l 

m=l 

m- 

M 

M 

M 

Yb^e  = 

m  m 

m=l 

m=\ 

m=\ 

0  ...  -a/2 

<x< 

M 

M 

M 

-  H<V2 


m=l 


m=\ 


m=\ 


Note  that  as  a  consequence  of  the  uniqueness  theorem,  it  is  only  necessary  to  satisfy  the 
magnetic  field  boundary  condition  in  the  region  |x|  <  t/2  since  the  tangential  electric  field 


boundary  condition  is  known  absolutely  (i.e.  =  0)  in  the  regions  -a/2  <  x  <  -t/2 


and  t/2<x<a/2  at  z  =  0,t/^. 

Enforcement  of  the  boundary  conditions  across  the  material  interfaces  in  the 
sample  holder  region  is  inherently  built  into  the  wave  transmission  analysis  outlined  in 
the  previous  chapter.  In  addition,  as  a  result  of  the  mode  orthogonality  that  exists  within 
the  sample  holder  region,  the  amplitude  coefficients  in  the  first  region  are  related  to  those 
in  layer  N  by 


1 _ 

_yij 

i=2  [ 

1 _ _ 

1 

3 

+ 
s 


^m.N 


^  m.N 


-^12 

Am 

^22 


m.N 


^  m,N 


(3.19) 


where  now  the  c.,b.  of  (2.35)  are  replaced  by  and 
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r"’'  =l  +  r’'' 

m  m 


(3.20) 


Z"  '  -Z 

m _ m 

Zlf  +  Zl 


The  index  m  expresses  that  the  A-parameters  will  be  different  for  eaeh  mode  that  is 
exeited  at  the  waveguide  to  sample  holder  junetions.  However,  sinee  the  height  and 
width  remain  eonstant  throughout  the  entire  sample  holder,  eaeh  mode  will  eouple 
independently  though  eaeh  layer,  as  mentioned  in  the  mode  orthogonality  diseussion 
above.  For  example,  if  the  TE30  mode  is  present  in  layer  one,  then  it  will  eouple 
eompletely  into  the  TE30  mode  in  layers  two,  three,  four,  ete.,  but  will  not  affeet  the  TEio 
or  TE50  modes  in  those  layers. 


3.3  Testing  with  Integral  Operators 

Before  implementing  the  last  step  of  the  modal-analysis  proeess,  it  is  benefieial  to 
review  the  results  of  the  previous  boundary  eondition  enforeement  seetion.  If  j ,  5^  j 
from  (3.19)  are  substituted  into  (3.15)-(3.16)  and  then  eaeh  remaining  eoeffieient  in 
(3.15)-(3.18)  is  divided  by  ,  there  will  be  4  linearly  independent  equations  with  4M 

unknown  ratios  (having  numerators  )•  However,  to  solve  for  these 

unknown  eoeffieient  ratios,  it  is  neeessary  to  have  AM  linearly  independent  equations 
with  AM  unknowns.  Similar  to  a  method  of  moments  teehnique  [10],  applieation  of 
properly  ehosen  inner  produet  testing  operators  will  generate  the  required  AM  equations. 
Consider  the  testing  operators 

i/2  a/2 

I  I  e„-{  ^dxdy  ,  n  =  l,2,...M  (3.21) 

-i/2  -a/2 
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I  ■{  ]dxdy  ,  n  =  \,2,..M 

-bjl  -tjl 


(3.22) 


where  the  uniqueness  theorem  is  again  used  to  limit  the  magnetie  field  intensity  operator 
to  -t/2  <  X  <  t/2  and  Z*  multiplies  the  inner  produet  to  express  the  faet  the  operator  is 
not  linked  to  a  speeifie  material  layer  in  the  sample  holder  region.  Applying  (3.21)  to 
(3. 15), (3. 17)  and  (3.22)  to  (3. 16), (3. 18)  produees 


a/2  M  V2  M  ^  I  72  M  ^  ^ 

f  e  ■  e.dx  +  V  [  e  ■  e  dx  =  V  \e- e'dx  +  V  I  e  -  e'dx  (3.23) 

V  \  n  1  ^  J  ^  J  ^  \  n  m  x  ■' 

0  m=l  0  m=\  0  ^=1  0 


tjl  M  72 

2,^ z;  f  r  •  kdx -Ta^ zi  f  r  •  hjx 

1  «  I  «  1  /  j  m  n  }  n  m 

0  “=1  0 


M  72  ^  ^  M  72  ^  ^ 

=  V  5^  .e''” Z*  f  h'  ■  h^’^dx -  V  Z*  f  •  h^^dx 

m,i  n  }  n  m  m,i  n  J  n  m 


(3.24) 


M  «/2  M  72  M  72 

j  =  j  •e;^Jx  +  £5„^  j  e„  -e^Jx 

m=l  0  m=l  0  m=l  0 


(3.25) 


M  72 


f  K  ■  TV  f  ■  hm’^dx “  V  s  ^  Z'  {  k'  ■  dx  (3 .26) 

wt  «  J  «  W2  m,i\  n  }  n  m  m,A«J«m  x  ■' 


m=\  0 


where  the  integration  over  y  eaneels  out  due  to  the  y-invarianee  of  the  problem  and 
evenness  in  x  was  also  used  to  limit  the  integration  to  only  half  the  waveguide  or  sample 
holder  width.  The  resulting  integrals  have  the  elosed-form  solutions 


C  =  \  e  -e  dx  =  \  A 

nm  J  -r 


—  ...  m  =  n 


021) 


0  ...  m^n 
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D  =  [  e  ■  dx  ■ 

nm  \  n  m 
0 


—  k'  =  k 

^  *  *  *  ^xm  ^xn 

(  rs  r 


(3.28) 


t/2  42 

F  =Z^  [¥ -h  dxM 

nm  n  ]  n  m  j  ^ 

"  (-ir¥-. 


...  t., =fc’ 


«  A:1  cosM 


Z  k^  —  k^^  ' ' ' 

■^W  ^x« 


(3.29) 


P  =r  ■h^’^dx  = 

nm  n  \  n  m 


...  m  =  n 


0  ...  m^n 


(3.30) 


Q  =Z' \h' -h^’^dx 

^nm  n  }  n  m 


...  m  =  n 


0  ...  m^n 


(3.31) 


where  the  details  are  loeated  in  Appendix  A  for  the  interested  reader. 

Substitution  of  from  (3.19)  into  (3 .23), (3. 24)  and  dividing  the  remaining 


eoeffieients  in  (3.23)-(3.26)  by  gives 


M  M  '  '  ~  ^  M  '  '  ~  ? 

c„,  =  £?■.  )  +  £'■-  )  (3.32) 


F  -'y  R  F  =yT  P  ( r  P  +  A'"e^^"‘‘‘' 


(3.33) 


IVI  IVl  IVl 

y  R  C  =yT  D  +y  r  D 

m  nm  m  nm  m  n. 


(3.34) 


y  R  F  =yTQ  -Yr  Q 

m  nm  m^nm  m^n, 


(3.35) 
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with  the  coefficient  ratios 


a 

R 


T.=- 


K..  =  ■ 


^  _  m 

m  + 


(3.36) 


One  final  rearrangement  leads  to  the  desired  result  of  AM  linearly  independent  equations 
with  AM  unknowns: 


VC  R  -Yu  T  -YV  r  =-C 

nm  m  nm  m  nm  m  h 

m=l  m=\  m=\ 

M  M  M 

ZF  R  +YW  T  -Yy  r  =F 

nm  m  nm  m  nm  m  n\ 

i=l  m=l  m=\ 

M  M  M 

Yd  t  +Y d  r  -Yc  t  =0 

nm  m  nm  m  nm  m 

m=\  m=\  m=\ 

M  M  M 

YQ  T  -YQ  r  -Y F  t  =0 

^nm  m  ^nm  m  nm  m 


m=l 

M 


n  =  \,2,...,M 


which  can  then  be  written  in  the  partitioned-matrix  form 


c 

-u 

-V 

0  ■ 

~R~ 

X, 

F 

w 

-Y 

0 

T 

X2 

0 

D 

D 

-c 

r 

0 

0 

Q 

-Q 

-F 

t 

0 

(3.37) 


(3.38) 


where  C,U,V,F,W,Y,D,Q  are  MxM  sub-matrices  and  are  the  column  vectors 


given  as 


Xi  =- 

,  X2  = 

^21 

1 

1 ..  _ 

(3.39) 
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3.4  Material  Parameter  Extraction. 


The  matrix  equation  of  (3.38)  is  readily  solved  for  the  coefficient  ratios  of  (3.36) 
by  Gaussian  elimination.  Since  the  waveguide  will  be  setup  for  dominant  TEio  mode 
operation,  all  the  higher-order  reflected  and  transmitted  modes  excited  by  the  reduced- 
width  sample  holder  will  evanesce  before  returning  to  the  network  analyzer.  That  is, 
while  the  higher-order  modes  will  influence  the  TEio  propagating  mode  amplitudes,  only 
the  reflected  and  transmitted  components  of  the  TEio  mode  will  arrive  at  the  detector 
elements.  Thus,  the  theoretical  coefficient  ratios  of  primary  interest  are 

R,=Sf^{(o,s,^)  ,  (3.40) 

To  extract  the  material  parameters  at  each  frequency,  an  iterative  optimization 
algorithm  is  used  to  find  the  complex  permittivity  and  permeability  that  best  match  the 
theoretically  calculated  to  the  experimentally  measured  S-parameter  as  described  by 
(3.1).  Eor  a  homogenous  material,  a  nonlinear  least-squares  algorithm  utilizing  all  four 
S-parameters  from  the  forward  and  reverse  traveling  waves  perhaps  offers  the  best 
method  to  filter  out  the  noise  inherent  in  the  system.  However,  a  complex  2-dimensional 
Newton's  root-searching  algorithm  can  be  used  to  solve  the  forward  and  reverse  problems 
separately  and  provide  insight  into  possible  inhomogeneities  or  gross  measurement 
errors. 

Note  that  the  material  parameter  extractions  presented  in  this  thesis  were  performed 
by  Matlab’s®  nonlinear  least-squares  algorithm.  As  a  precaution  to  the  reader,  this 
algorithm  requires  all  variables  to  be  real.  Thus,  it  is  necessary  to  break  the  permittivity 
and  permeability  into  their  respective  real  and  imaginary  parts. 
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3.5  Summary 


A  three  step  modal  analysis  teehnique  was  developed  to  aeeommodate  the  higher- 
order  modes  exeited  by  the  stepped  junetion  between  the  waveguide  feeds  and  the  multi¬ 
layered  redueed-width  sample  holder  region.  The  problem  was  formulated  into  a 
partitioned-matrix  form  to  solve  for  the  theoretieal  S-parameters  for  a  guess  at  the 
eonstitutive  material  parameters.  Finally,  two  optimization  algorithms  for  extraeting  the 
test  sample's  material  parameters  from  the  experimental  S-parameters  were  suggested. 
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4.  Results  and  Error  Analysis 


To  verify  the  modal-analysis  teehnique  developed  in  Chapter  3,  several  materials 
were  measured  in  the  S-band  frequency  range,  2.6-3.95  GHz  .  While  the  technique  is 
intended  for  low-frequency  applications  with  large  test  samples,  the  smaller  S-band 
waveguide  samples  were  adequate,  and  much  less  expensive,  for  a  proof-of-concept 
demonstration.  In  this  chapter,  experimental  results  will  be  shown  for  acrylic,  magnetic 
radar  absorbing  material  (MRAM),  and  electric  radar  absorbing  material  (ERAM).  In 
addition,  a  brief  study  of  the  number  of  higher-order  modes  required  for  material 
parameter  convergence  for  different  sample  holder  aperture  reductions  will  be  presented. 

When  conducting  any  experimental  procedure,  there  are  many  possible  error 
sources.  Assuming  a  good  calibration,  properly  machined  sample  holders  with  precision 
alignment  pins,  and  quality  test  samples  with  uniform  thickness  and  material 
composition,  user  error  can  still  lead  to  inaccuracies  in  the  constitutive  parameter 
extractions  when  performing  these  waveguide  measurements.  A  differential  error 
analysis  will  be  presented  to  develop  error  bounds  related  to  test  sample  thickness  and 
placement  measurement  errors. 

4.1  Experimental  Setup 

The  experimental  setup  is  as  shown  in  Figure  4-1  with  a  WR284  waveguide 
(where  WR  stands  for  “waveguide  rectangular”  and  284  refers  to  the  inner  waveguide 
width  a  =  2.84  inches  )  connected  to  a  HP85  IOC  Vector  Network  Analyzer  via  coaxial 
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cables.  The  redueed-aperture  sample  holder  is  mounted  to  the  waveguide  with  preeision 
alignment  pins  to  minimize  the  diseontinuities  aeross  the  interfaees.  Also,  great  eare  was 
taken  to  elean  all  eabling  junctions  with  isopropyl  alcohol  to  ensure  proper  eleetrieal 
eontaets.  Finally,  a  full  two-port  thru,  reflect,  line  (TRL)  ealibration  proeedure  was  used 
to  set  the  ealibration  planes  immediately  before  and  immediately  after  the  sample  holder 
(at  z  =  0,  with  referenee  to  Figure  3-1)  [9]. 

The  three  sample  holders  used  to  verify  the  modal-analysis  teehnique  are  shown 
in  Figure  4-2.  The  sample  holder  on  the  left  has  the  same  dimensions  as  the  waveguide 
feeds  and  was  used  to  set  the  baselines,  sinee  only  the  dominant  TEio  mode  was  required 
for  material  parameter  extraetion.  The  sample  holders  in  the  middle  and  on  the  right  of 
Figure  4-2  have  redueed  width,  but  no  reduetion  in  height.  It  will  be  shown  that 
eomparable  results  ean  be  obtained  using  the  redueed-width  sample  holders  by 
aeeounting  for  the  higher-order  modes  that  are  exeited  at  the  waveguide  feed  to  sample 
holder  junetions. 

Even  before  measuring  any  test  samples  there  were  a  few  additional  steps  taken  to 
validate  the  experimental  setup.  Eirst,  after  performing  the  ealibration,  the  shorting 
standard  (brass  plate)  used  in  the  reflect  portion  of  the  TRL  ealibration  proeedure  was 
measured  again  to  ensure  that  »  -1  -yO  and  »  0  -  yO  .  Also,  an  extraetion  was 

performed  on  the  S-parameters  measured  from  the  empty  sample  holders  to  ensure  a  free 
spaee  baseline  for  the  system. 
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Figure  4-1.  Experimental  setup  showing  a  WR284  waveguide  connected  to  an  HP8510C 
Vector  Network  Analyzer  with  coaxial  cables.  The  reduced-aperture  sample  holder  is 
mounted  between  the  waveguide  feeds  with  precision  alignment  pins. 


Figure  4-2.  S-band  waveguide  sample  holders.  On  the  left  is  the  full-aperture  sample 
holder  having  the  same  dimensions  as  the  waveguide  feeds,  while  the  ones  in  the  middle 
and  on  the  right  have  a  reduced  width,  but  no  reduction  in  height. 
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4.2  Error  Analysis 


Before  proceeding  with  the  differential  error  analysis  development,  it  is  necessary 
to  first  review  the  big  picture  goal  of  the  modal-analysis  technique  derivation  of  the 
previous  chapter.  The  complex  relative  permittivity  and  permeability  of  the  test  sample 
will  be  determined  using  an  iterative  nonlinear  least-squares  optimization  algorithm  that 
seeks  to  match  the  theoretically  calculated  to  the  experimentally  measured  S-parameters: 


minJ^[S'f^(®,f,//) 


(3.1) 


where  i  =  \  1,21,22,12  .  Thus,  uncertainties  in  the  experimental  S-parameters  will  clearly 
lead  to  different  s,ili  values  that  satisfy  (3.1).  In  addition,  given  “perfect”  experimental 
S-parameters  with  no  added  random  noise  or  drift,  there  are  obviously  many  other 
possible  error  contributions  that  will  lead  to  extraction  inaccuracies.  Uncertainties  in 
sample  holder  dimensions  (width,  height,  length  -  t,  b,  dj^ ),  alignment  errors  causing 

excitation  of  undesired  higher-order  modes,  and  absorption  due  to  the  imperfectly 
conducting  sample  holder  walls  are  just  a  few  possible  error  sources.  However,  carefully 
machined  sample  holders  with  precision  alignment  pins  lead  to  very  good  empty-aperture 
baselines,  and  experience  has  shown  that  user  error  in  measuring  test  sample  thickness 
and  placement  within  the  sample  holder  are  two  of  the  most  sensitive  and  difficult  to 
control  parameters. 

To  account  for  the  test  sample  thickness  and  placement  measurement  errors,  a 
differential  error  analysis  was  performed.  With  the  estimated  permittivity  written  as  a 
function  of  thickness  I  and  placement  p  , 
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£  =  s'{i,p)-j£"{i,p) 


(4.1) 


the  Taylor  series  expansion  of  the  real  part  around  the  point  (f  is  given  as 

£'{£^+A£,p^+Ap)  =  £'{£^,p^)+^^  + 

o£  op 


(4.2) 


where  A£  is  the  uneertainty  in  the  measured  thiekness  £^ ,  and  Ap  is  the  uneertainty  in 
plaeement  at  the  loeation  p^ .  For  relatively  small  A£  and  Ap ,  the  higher-order  terms  of 

the  Taylor  series  expansion  ean  be  negleeted  and  the  ehange  in  the  real  permittivity  ean 
be  represented  as 


=  £'(«,  +  A<,p„  +  Ap)-£'(«„,;>„)*  A<  +  Ap  =  (4.3) 


Next,  the  partial  derivatives  ean  be  approximated  (again  assuming  small  A£  and  Ap)  by 


A£ 


dp 


Ap 


In  general,  A£  and  Ap  ean  be  positive  or  negative  values.  For  small  relative 


perturbations,  the  respeetive  errors  and  A^'  will  be  symmetrie  about  the 


unperturbed  £'[£^,p^) .  However,  for  larger  perturbations,  the  symmetry  will  start  to 

break  down  and  the  errors  generated  by  the  positive  and  negative  values  must  be 
eomputed  separately  to  determine  the  worst-ease  eontributions. 

It  should  also  be  elear  that  independent  material  parameter  extraetions  for  the 
perturbations  in  the  measured  thiekness  and  plaeement  of  the  test  sample  ean  be 
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conducted  and  then  eombined  to  determine  the  overall  error  bounds  related  to  those 
speeifie  uneertainties.  In  eombining  the  results,  the  triangle  inequality  is  used  to  form  an 
approximation  of  the  worst-ease  errors  in  the  real  permittivity 

\^<ppro}\  =  +  (4-6) 

The  triangle  inequality  is  generally  applied  to  veetors  with  the  geometrie  interpretation 
that  the  sum  of  the  lengths  of  any  two  sides  of  a  triangle  is  always  greater  than  the  length 
of  the  remaining  side.  In  this  ease,  where  lS.s[,/S.s'p  are  sealar  quantities,  the  inequality  is 

simply  used  to  prevent  eaneellation  between  positive  and  negative  error  eontributions. 

It  must  be  noted  that  often  independent  error  eontributions  are  added  in 

quadrature  [4]  (i.e.,  |A^'|  =  yj(Asff+~(As^  ),  but  that  assumes  normally  distributed 

errors  where  the  sum  of  all  deviations  around  the  nominally  measured  values  are 
negligible.  In  addition,  the  least-squares  optimization  algorithm  tends  to  further 
exaeerbate  the  problem,  thus  favoring  the  worst-ease  approaeh.  Finally,  the  error 
analyses  for  the  imaginary  permittivity  and  real  and  imaginary  permeability  follow  in  a 
similar  manner. 

4,3  Acrylic 

The  aerylie  measured  for  this  thesis  was  a  elear  plexiglass-type  material  that 
exhibited  almost  pure  dieleetrie  eharaeteristies  with  very  little  eleetrie  loss  and  no 
magnetie  properties.  In  addition,  due  to  its  rigid  strueture,  aerylie  is  easily  maehined  into 
tightly  fitting  waveguide  samples  with  almost  uniform  thiekness.  Beeause  of  these 
favorable  qualities,  aerylie  was  an  ideal  material  to  start  the  modal-analysis  verifieation. 


4-6 


Having  a  priori  knowledge  of  the  dieleetrie  nature  of  the  test  sample,  a 
permittivity  only  extraetion  using  all  four  S-parameters  was  eondueted  first  on  a  full-size 
test  sample  to  set  the  baseline  relative  dieleetrie  eonstant  at  about  2.65  aeross  the  S-band 
frequeney  range.  Note  that  the  test  sample  was  mounted  flush  toward  the  front  surfaee  of 
the  full-aperture  sample  holder  (faeing  terminal  one  of  the  network  analyzer)  to  minimize 
the  diffieulty  in  measuring  sample  plaeement. 

A  seeond  aerylie  test  sample  was  measured  in  the  2.34-ineh  redueed-width 
sample  holder.  Figure  4-3  shows  that  about  10  modes  were  required  to  eonverge  to  the 
full-aperture  baseline.  That  is,  the  energies  spread  into  higher-order  modes  must  be 
aeeounted  for  to  obtain  aeeurate  results.  Note  that  error  bounds  related  to  uneertainty  in 
test  sample  thiekness  or  plaeement  were  not  eonsidered  sinee  the  results  were  virtually 
identieal  to  the  full-aperture  baseline. 

Finally,  a  third  aerylie  test  sample  was  measured  in  the  1 .34-ineh  redueed-width 
sample  holder.  With  a  larger  disruption  in  the  waveguide  wall  surfaee  eurrents,  more 
modes  were  required  for  eonvergenee  to  the  full-aperture  baseline.  As  shown  in 
Figure  4-4,  about  20  modes  were  required  to  obtain  aeeurate  results  in  the  more  tightly 
choked  sample  holder.  Note  that  the  imaginary  permittivity  graphs  were  not  explicitly 
shown  since  they  exhibited  very  low  loss  and  added  nothing  additional  to  the  discussion. 
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a  =  2.84  in  t=  2.34  in 
L  =  118.8  +  0  miis  AP  =  0  miis 
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Figure  4-3.  Real  permittivity  for  acrylic  test  sample  mounted  flush  toward  port  one  in  the 
2.34-inch  sample  holder.  About  10  modes  are  required  for  convergence  to  the  full- 
aperture  baseline.  (Permittivity  only  extraction  using  least-squares  algorithm  on  all  four 
S-parameters.) 


a  =  2.84  in  t=  1.34  in 
L=  119  +  0  miis  AP  =  0  miis 
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Figure  4-4.  Real  permittivity  for  acrylic  test  sample  mounted  flush  toward  port  one  in  the 
1.34-inch  sample  holder.  About  20  modes  are  required  for  convergence  to  the  full- 
aperture  baseline.  (Permittivity  only  extraction  using  least-squares  algorithm  on  all  four 
S-parameters.) 
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4.4  MRAM 


Moving  on  to  a  more  complicated  extraction,  three  properly  sized  MRAM  test 
samples  were  measured  employing  the  same  basic  procedures  used  for  the  acrylic.  In  all 
cases,  the  samples  were  mounted  flush  to  the  front  surfaces  of  the  sample  holders  to 
minimize  placement  uncertainty.  In  addition,  when  using  the  reduced-width  sample 
holders,  the  number  of  modes  was  increased  until  convergence  was  achieved.  However, 
since  the  MRAM  material  was  expected  to  exhibit  magnetic  properties,  the  nonlinear 
least-squares  algorithm  was  setup  for  both  a  complex  relative  permittivity  and 
permeability  extraction. 

The  results  for  the  MRAM  measured  in  the  2.34-inch  reduced-width  sample 
holder  are  compared  to  the  full-aperture  baselines  in  Figure  4-5  and  Figure  4-6.  The 
permeability  is  seen  to  be  nearly  identical  with  the  full-aperture  baseline,  and  the  minor 
difference  in  the  real  permittivity  is  readily  explained  by  a  2-mil  uncertainty  in  thickness 
and  a  1-mil  uncertainty  in  placement.  The  MRAM  material  used  was  manufactured  with 
an  adhesive  on  one  side,  which,  when  removed  with  sandpaper,  led  to  slight  fluctuations 
in  thickness  across  the  sample. 

Similarly,  Figure  4-7  and  Figure  4-8  show  the  results  for  the  MRAM 
measurement  in  the  1.34-inch  sample  holder.  Again,  the  permittivity  differences  are 
mostly  explained  by  the  thickness  and  placement  uncertainties.  However,  the  error  bars 
do  not  completely  account  for  the  1.45%  maximum  discrepancy  in  the  real  permeability 
relative  to  its  full-aperture  baseline.  This  small  discrepancy  is  mainly  attributed  to  the 
grooves  left  by  the  sandpaper,  but  other  possible  errors  due  to  uncertainty  in  the 
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experimental  S-parameters  and  sample  holder  dimensions  and  alignment  likely 
contributed  to  the  deviation  from  the  full-aperture  real  permeability  baseline.  Better 
results  would  most  likely  be  attained  with  a  test  sample  that  was  manufactured  without  an 
adhesive,  or  if  a  less  destructive  method  was  used  to  remove  the  adhesive.  In  addition, 
each  test  sample  was  cut  separately  from  a  large  sheet.  Improvements  might  be  made  if 
the  full-aperture  test  sample  were  measured,  and  then  cut  down  for  a  measurement  in  the 
2.34-inch  sample  holder,  and  then  cut  down  again  for  a  measurement  in  the  1.34-inch 
sample  holder.  By  using  portions  of  the  original  full-aperture  test  sample,  errors  due  to 
fluctuations  in  material  composition  across  the  large  sheets  would  be  reduced. 
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a  =  2,84  in  t  =  2.34  in  l\/l  =  15  modes 
L  =  187  +  2  mils  AP  =  1  mils 


a  =  2.84  in  t  =  2.34  in  M  =  1 5  modes 
L  =  1 87  ±  2  mils  AP  =  1  mils 


.  Reduced  Aperture 

-  Full  Aperture 


Figure  4-5.  Permittivity  of  MRAM  test  sample  mounted  flush  toward  port  one  in  the 
2.34-inch  sample  holder.  Small  difference  in  real  part  easily  explained  by  uncertainty  in 
measured  thickness  and  placement. 


a  =  2.84  in  t  =  2.34  in  M  =  1 5  modes 
L  =  187  +  2  mils  AP  =  1  mils 


a  =  2.84  in  t  =  2.34  in  M  =  15  modes 
L  =  187  +  2  mils  AP  =  1  mils 

5 1 - . - . - ^ - . - . - ^ 


4 


3 


1 

qI - 1 - r - ^ - r - 1 - ^ 

2.6  2.8  3  3.2  3.4  3.6  3.8 

Frequency  (GHz) 


.  Reduced  Aperture 

-  Full  Aperture 


Figure  4-6.  Permeability  of  MRAM  test  sample  mounted  flush  toward  port  one  in  the 
2.34-inch  sample  holder. 
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a  =  2,84  in  t  =  1 .34  in  M  =  20  modes 
L  =  1 98  +  2  mils  AP  =  1  mils 


Frequency  (GHz) 


a  =  2.84  in  t  =  1 .34  in  M  =  20  modes 


L  =  1 98  ±  2  mils  AP  =  1  mils 


Frequency  (GHz) 


.  Reduced  Aperture 

-  Full  Aperture 


Figure  4-7.  Permittivity  of  MRAM  test  sample  mounted  flush  toward  port  one  in  the 
1.34-inch  sample  holder.  Thickness  and  placement  uncertainty  explain  most  of  the 
difference  relative  to  the  full-aperture  baseline. 


a  =  2.84  in  t  =  1 .34  in  M  =  20  modes 
L  =  198±  2  mils  AP  =  1  mils 


a  =  2.84  in  t  =  1 .34  in  M  =  20  modes 


L  =  198±  2  mils  AP  =  1  mils 


.  Reduced  Aperture 

-  Full  Aperture 


Figure  4-8.  Permeability  of  MRAM  test  sample  mounted  flush  toward  port  one  in  the 
1.34-inch  sample  holder.  Maximum  real  permeability  difference  is  1.45%  of  full- 
aperture  baseline  without  considering  the  thickness  and  placement  uncertainties. 
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4.5  ERAM 


The  final  material  measured  for  this  thesis  was  the  ERAM,  whieh  consisted  of  a 
carbon-based  foam  from  which  it  proved  very  difficult  to  cut  precision  test  samples. 
Loose-fitting  test  samples  led  to  energy  leaking  through  the  air  gaps  between  the  test 
sample  and  the  sample  holder  walls.  On  the  other  hand,  test  samples  that  were  cut  too 
large  caused  compression  of  the  foam  also  leading  to  inaccuracies.  In  addition,  it  was 
very  difficult  to  accurately  measure  the  thickness  of  the  test  samples  since  the  foam  was 
easily  compressed  by  the  inch  micrometer. 

Three  separate  attempts  were  made  to  measure  ERAM.  The  results  shown  in 
Eigure  4-9  and  Eigure  4-10  are  representative  of  all  three  attempts.  The  error  bars  related 
to  a  5-mil  uncertainty  in  measured  thickness  and  a  1-mil  uncertainty  in  placement  do  not 
fully  explain  the  discrepancies  between  the  reduced-width  test  samples  and  the  full- 
aperture  baseline.  Similar  error  bars  added  to  the  full-aperture  baseline  would  still  not 
completely  account  for  the  differences.  Without  taking  the  error  bars  into  account,  the 
real/imaginary  permittivity  of  the  test  sample  measured  in  the  2.34-inch  sample  holder  is 
within  6.55/14.22%  of  the  full-aperture  baselines.  Similarly,  the  real/imaginary 
permittivity  of  the  ERAM  measured  in  the  1 .34-inch  sample  holder  is  within 
1 1.01/22.71%  of  the  full-aperture  baselines. 

When  extraction  results  exhibit  large  inconsistencies,  it  is  often  informative  to 
analyze  the  raw  experimental  S-parameters.  Table  4-1  displays  the  measured 
S-parameter  ranges  from  all  three  ERAM  test  samples.  Note  that  due  to  the  large  loss 
tangent  of  the  ERAM  material,  the  transmission  is  quite  low,  especially  for  the  1 .34-inch 
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test  sample.  With  only  4.6%  of  the  power  transmitted  through  the  test  sample,  the  noise 
floor  of  the  network  analyzer  beeomes  a  coneern.  An  extraetion  based  on  only  the 
refleetion  eoeffieients  ( 5*“^ ,  5*2“^ )  was  attempted,  but  the  results  were  not  improved. 

In  summary,  the  ineonsisteneies  in  the  ERAM  extraetion  results  are  attributed  to 
three  main  error  sourees.  First,  the  ineonsisteneies  highlight  the  need  for  very  high 
quality  waveguide  test  samples  that  fit  into  the  sample  holders  without  air  gaps  or 
eompression.  Seeond,  the  ERAM  foam  is  even  more  suseeptible  than  the  MRAM  to 
material  fluetuations  (i.e.,  inhomogeneities)  aeross  the  sheet,  and  better  results  would 
most  likely  be  aehieved  by  eutting  the  2.34  and  1.34-meh  test  samples  from  the  full- 
aperture  pieee  as  deseribed  in  the  previous  seetion.  And  third,  the  highly  loaded  radar 
absorbing  material  eauses  a  low  transmitted  signal-to-noise  ratio,  espeeially  in  the  tightly 
ehoked  1.34-meh  sample  holder. 


Table  4-1 .  Experimental  S-parameter  ranges  for  ERAM  (in  linear  power  format). 
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a  =  2,84  in  t  =  2.34  in  l\/l  =  15  modes 


L  =  240  +  5  mils  AP  =  1  mils 


Frequency  (GHz) 


a  =  2.84  in  t  =  2.34  in  M  =  15  modes 
L  =  240  ±  5  mils  AP  =  1  mils 


Frequency  (GHz) 


.  Reduced  Aperture 

-  Full  Aperture 


Figure  4-9.  Permittivity  of  ERAM  test  sample  mounted  flush  toward  port  two  in  the 
2.34-inch  sample  holder.  Maximum  difference  of  real/imaginary  permittivity  is 
6.55/14.02%  relative  to  their  full-aperture  baselines,  without  considering  thickness  and 
placement  uncertainties.  (Permittivity  only  extraction  using  all  four  S-parameters.) 


a  =  2.84  in  t  =  1 .34  in  M  =  20  modes 


L  =  240  ±  5  mils  AP  =  1  mils 


Frequency  (GHz) 


a  =  2.84  in  t  =  1 .34  in  M  =  20  modes 
L  =  240  ±  5  mils  AP  =  1  mils 


Frequency  (GHz) 


.  Reduced  Aperture 

-  Full  Aperture 


Figure  4-10.  Permittivity  of  ERAM  test  sample  mounted  flush  toward  port  two  in  the 

I. 34-inch  sample  holder.  Maximum  difference  of  real/imaginary  permittivity  is 

I I. 01/22.71%  relative  to  their  full-aperture  baselines,  without  considering  thickness  and 
placement  uncertainties.  (Permittivity  only  extraction  using  all  four  S-parameters.) 
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4.6  Summary 


Acrylic,  MRAM,  and  ERAM  test  samples  mounted  in  reduced-width  S-band 
sample  holders  were  used  to  verify  the  modal  analysis  teehnique.  The  aerylie  and 
MRAM  results  were  exeellent,  and  the  ERAM  results  were  reasonable  eonsidering  the 
difficulties  in  working  with  the  high-loss  earbon-based  foam.  It  was  shown  that  the 
number  of  higher-order  modes  required  for  convergence  is  directly  related  to  the  degree 
of  width  reduetion  to  the  sample  holder.  Einally,  a  differential  error  analysis  was 
presented  to  form  error  bounds  related  to  human  errors  in  measurement  of  test  sample 
thickness  and  placement. 
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5.  Conclusion  and  Recommendations 


A  stepped-waveguide  technique  was  developed  to  facilitate  permittivity  and 
permeability  measurement  of  smaller  test  samples  in  larger  rectangular  waveguide 
measurement  systems.  Specifically,  a  mode-matching  technique  was  implemented  to 
accommodate  higher-order  modes  excited  by  the  stepped  junction  between  the  waveguide 
feeds  of  the  measurement  system  and  the  reduced-width  sample  holder.  The  sample 
permittivity  and  permeability  were  estimated  by  using  a  nonlinear  least-squares  algorithm 
to  compare  the  theoretical  S-parameters  obtained  from  the  modal  analysis  and  the 
experimentally-measured  S-parameters  obtained  from  the  network  analyzer. 

The  reduced-aperture  material  characterization  technique  was  experimentally 
verified  through  comparisons  of  stepped-waveguide  measurements  of  acrylic,  MRAM, 
and  ERAM  with  traditional  full-aperture,  single-mode  measurements.  The  comparisons 
demonstrated  the  validity  and  accuracy  of  the  technique  and  also  revealed  that  smaller 
apertures  require  a  greater  number  of  modes  to  achieve  convergence,  as  anticipated.  A 
differential  error  analysis  conducted  on  test  sample  thickness  and  placement  uncertainties 
explained  most  of  the  discrepancies  from  the  full-aperture  baselines.  Finally,  the  high- 
loss  ERAM  demonstrated  the  limitations  of  the  technique  when  there  is  poor  transmission 
coupling  combined  with  loose-fitting  or  compressed  test  samples. 
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5.1  Future  Research 


While  experimental  data  from  low-frequency  waveguide  measurements  would 
perhaps  offer  insight  into  the  best  direction  for  future  investigation,  there  are  a  few 
general  improvements  that  could  be  made  to  the  current  reduced-width  modal-analysis 
technique.  In  addition,  several  new  modal-analysis  techniques  warrant  investigation. 

5.1.1  Possible  Improvements  to  Current  Technique. 

First,  the  differential  error  analysis  could  be  expanded  to  include  the  effects  of 
uncertainties  in  more  variables,  such  as  the  experimental  S-parameters  and  sample  holder 
dimensions.  However,  each  additional  variable  uncertainty  requires  a  separate  material 
extraction  to  obtain  its  related  error  bounds,  which  are  then  added  to  the  effects  of  the 
other  variables  using  the  triangle  inequality. 

Next,  from  a  numerical  standpoint,  the  most  time-consuming  part  of 
implementing  the  technique  is  solving  matrix  (3.38)  for  the  theoretical  S-parameters 
(3.40).  Investigation  of  a  more  efficient  algorithm  than  the  simple  Gaussian  elimination 
method  used  for  this  thesis  would  expedite  the  extraction  process.  As  a  reference,  using  a 
2-GHz  computer  processor,  a  typical  20-mode  least-squares  extraction  (involving  all 
four  S-parameters)  of  201  frequency  points  for  a  1. 34-inch  wide  test  sample  takes  about 
2-3  minutes.  Then,  each  variable  uncertainty  requires  about  the  same  extraction  duration. 
Thus,  while  each  extraction  is  by  no  means  an  all  day  process,  analyzing  several 
measurements  can  develop  into  a  tedious  endeavor. 

Finally,  an  attenuation  correction  factor  for  the  imperfectly  conducting  sample 
holder  walls  could  be  explored  using  a  perturbation  approach.  Another  correction  factor 


5-2 


might  also  be  found  for  air  gaps  between  the  test  sample  and  the  sample  holder  walls  by  a 
similar  method. 

5,1,2  Possible  New  Modal-Analysis  Techniques, 

A  dual  aperture  reduetion  (in  width  and  height)  would  allow  for  an  additional 
degree-of-freedom  in  the  maehining  of  test  samples.  The  aceounting  proeess  for  the 
possible  excited  mode  shapes  across  the  new  junction  geometry  would  have  to  be 
reinvestigated. 

Also,  a  modal-analysis  might  be  implemented  for  thin  resistive  sheets  (R-cards) 
clamped  directly  between  waveguide-feed  flanges  without  the  use  of  a  sample  holder. 
Copper  tape  could  be  applied  around  the  flanges  to  prevent  energy  from  escaping  from 
the  R-card  to  air  interfaces. 

Due  to  its  higher  resolution,  a  reduced-aperture  probe  might  aid  in  detecting  sub¬ 
surface  defects  in  low-observable  (LO)  coatings  for  low-frequency  radar  threats.  Instead 
of  determining  the  material  parameters,  the  probe  would  scan  the  surface  while  a  real¬ 
time  oscilloscope  trace  was  monitored  for  sharp  spikes. 

On  the  other  end  of  the  spectrum,  waveguide  test  sample  sizes  become 
exceedingly  small  for  high-frequency  applications.  When  dealing  with  such  small  test 
samples,  machining  inaccuracies  and  sample  uniformity  are  important  concerns.  A 
modal-analysis  technique  developed  for  a  dual  aperture  expansion  would  help  alleviate 
these  concerns  by  allowing  a  more  substantial  test  sample  to  be  measured  and  still  extract 
the  material  parameters  at  the  high  frequencies  of  interest. 
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Appendix  A.  Proof  of  Matrix  Relations 


This  appendix  contains  the  proofs  to  the  integral  solutions  given  in  (3.27)-(3.31). 
The  following  trigonometric  identities  are  used: 

2cosxcos  j  =  cos(x  + j)  +  cos(x- j)  (A.l) 


sin(x±  j)  =  sinxcos  j±cosxsin  j 


Also,  as  a  reminder,  in  the  waveguide  regions. 


e  =  vcos  A:  x  ,  h  = 
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and  for  the  layer  in  the  sample  holder  region. 
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